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ABSTRACT 

Aims. PSR J0205+6449 is a young rotation-powered pulsar in SNR 3C 58. It is one of only three young (< 10, 000 year old) pulsars 
which are so far detected in the radio and the classical X-ray bands, as well as at hard X-rays above 20 keV and at high-energy (> 100 
MeV) y-rays. The other two young pulsars are the Crab and PSR B 1509-58. Our aim is to derive the timing and spectral characteristics 
of PSR J0205+6449 over the broad X-ray band from ~ 0.5 to ~ 270 keV. 

Methods. We used all publicly available RXTE observations of PSR J0205+6449 to first generate accurate ephemerides over the 
period September 30, 2000 - March 18, 2006. Next, phase-folding procedures yielded pulse profiles using data from RXTE PC A 
and HEXTE, and XMM-Newton EPIC PN. All profiles have been phase aligned with a radio profile derived from the Jodrell Bank 
Observatory data, and the time-averaged timing and spectral characteristics of the pulsed X-ray emission have been derived. 
Results. While our timing solutions are consistent with earlier results, our work shows sharper structures in the PCA X-ray profile. 
The X-ray pulse profile consists of two sharp pulses, separated in phase by 0.488 ± 0.002, which can be described with 2 asymmetric 
Lorentzians, each with the rising wing steeper than the trailing wing, and full-width-half-maximum 1.41 ± 0.05 ms and 2.35 ± 0.22 
ms, respectively. We find an indication for a flux increase by a factor ~ 2, about 3.5cr above the time-averaged value, for the second, 
weaker pulse during a two-week interval, while its pulse shape did not change. The spectrum of the pulsed X-ray emission is of 
non-thermal origin, exhibiting a power-law shape with photon index r = 1.03 ± 0.02 over the energy band ~ 0.5 to ~ 270 keV. In the 
energy band covered with the PCA (~ 3 - 30 keV) the spectra of the two pulses have the same photon index, namely, 1.04 ± 0.03 and 
1 . 10 ± 0.08, respectively. Comparisons of the detailed timing and spectral characteristics of PSR J0205+6449 in the radio, hard X-ray 
and gamma-ray bands with those of the Crab pulsar, PSR B 1509-58 and the middle-aged Vela pulsar do reveal more differences than 
similarities. 

Key words. Stars: neutron - pulsars: individual PSR J0205+6449, PSR B1509-58, Crab, Vela - X-rays: general - Gamma rays: 
observations - Radiation mechanisms: non-thermal 



1. Introduction 

PSR J0205 +6449 is a young rotation-powered pulsar of which 
the pulsations were first discovered in X-rays in a 2002 
Chandra X-r a y Ob servatory (CXO) observation, reported by 
iMurrav et al.l d2002h together with a confirmation in an 
analysis of archival Rossi X-Ray Timing Explorer (RXTE) 
data. Subsequently, the weak radio signal was detected by 
ICamilo et al.l 0002). PSR J0205 +6449 is a young, 65-ms pul- 
sar located in the center of supernova remnant/pulsar wind neb- 
ula (PWN) 3C 58. It is one of the most energetic pulsars in the 
Galaxy with a spin-down luminosity E ~ 2.7 x 10 37 erg s _1 , 
and characteristic age r ~ 5.4 kyr. This characteristic age, es- 
timated with the values of the period and period derivative, puts 
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in doubt the possible association with 3C 58, which coincides 
positionally with the histori cal 828 yr old supernova SN1191 
(Stephens on & Greed. I20 02). However, an age of several thou- 
sand years for 3C 58, closer to the characteristic age of the pul- 
sar, can be derived from the velocities of the radio expansion of 
the P WN ( Biet enhoTzl I2006) and of optical knots dFesen et all 
I2008h . 

Recently, Livingst one et alJ ((2009) presented for 
PSR J0205 +6449 phase-coherent timing analyses using X- 
ray data from the Proportional Counter Array (PCA; 2-60 keV) 
aboard RXTE and radio data from the Jodrell Bank Observatory 
and the Green Bank Telescope (GBT), spanning together 6.4 
yrs. This work revealed timing noise and two spin-up glitches. 
Furthermore, they presented detailed characteristics of the X-ray 
profile, which was detected up to ~40 keV. Their X-ray profile 
template consisted of two Gaussian-shaped pulses, a narrow 
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(full-width-half-maximum (FWHM) ~ 1.6 ms), more intense 
pulse and a broader (FWHM -3.8 ms) weak pulse separated 
0.505 in phase, the single radio pulse leading the main X-ray 
pulse by = 0.10 + 0.01. Earlier results from an analysis of 
part o f the RXTE and GBT data were reported bv lRansom et alj 
d2004l) . These authors also presented spectral fits over the 
energy band 3-16 keV for both pulses: the best fit power-law 
photon indices were hard, namely F = 0.84^j ( ^ for the main 
pulse and F = 1 .0^ 3 for the second (weaker) pulse. 

Finally, high-energy y-ray pulsations (> 0.1 GeV) from 
PSR J0205 +6449 were discovered with the Large Area 
Telescope (LAT) ab oard the Fermi Gamma-ray Space Telescope 
dAbdo et al.L |2009a), folding the y-ray arrival times with the ra- 
dio rotational ephemeris from, again, the GBT and Jodrell Bank. 
The y-ray light curve for energies > 0.1 GeV shows also two 
peaks with intensities differing by a factor ~ 2, aligned with the 
X-ray peaks, However, the main X-ray pulse coincides in phase 
with the weakest y-ray pulse which has the softest spectrum of 
the two at high-energy y-rays. The total pulsed y-ray spectrum 
exhibits a simple power-law shape with index F ~ 2. 1 and expo- 
nential cutoff at ~ 3.0 GeV. 

PSR J0205+6449 is now one of only three young (< 10, 000 
year old) pulsars which are detected in the classical X-ray band 
and at hard X-rays above 20 keV, as well as at high-energy (> 0.1 
GeV) y-ray energies, the others being the Crab pulsar and PSR 
B1509-58 (PSR J1513-5908). The Crab pulsar has been studied 
over the total high-energy band already in great detail (see for a 
cohere nt high-energy pic ture from soft X-rays up to high-energy 
y-rays iKuiper et all 1200 with eve n a detection of pulsed y- 
rays above 25 GeV dAliu et alll2008l) . 

The detection of pulsed emission above 100 MeV from PSR 
B 1509-58 had to wait f or the new generation o f currently oper- 
ational y-ray telescopes (Pell izzoni et alll2009l) . However, these 
three young pulsars have very different timing and spectral char- 
acteristics. This makes it particularly interesting to determine the 
timing and spectral characteristics of PSR J0205 +6449 in more 
detail over the high-energy band of the electro-magnetic spec- 
trum in order to compare these with those of Crab and PSR 
B 1509-58 and for confrontation with theoretical predictions. In 
this work our aim is to extend the coverage in the hard X-ray 
band to higher energies, exploiting the data of the High Energy 
X-ray Timing Experiment (HEXTE; 1 5-250 keV) aboard RXTE, 
and to extend the energy window to lower energies by analysing 
data from XMM-Newton. We will present the results from our 
timing study exploiting only the multi-year PCA/RXTE moni- 
to ring data, which we perfo rmed in parallel to the work reported 
bv lLivingstone et all (120091) . Our timing solutions are consistent 
with those of the latter authors, but our work revealed sharper 
structures in the PC A X-ray pulse profile. Furthermore, we de- 
rive the spectral characteristics over the total X-ray band. In the 
discussion we compare our findings with the characteristics of 
PSR J0205 +6449 reported in the radio band and at high-energy 
y-rays, as well as with the timing and spectral characteristics of 
the Crab pulsar, PSR B1509-58 and the middle-aged Vela pulsar. 

2. Instruments and observations 

2.1. RXTE 

In this study extensive use is made of data from monitoring 
observations of PSR J0205 +6449 with the two non-imaging X- 
ray instruments aboard RXTE, the Proportional Counter Array 
(PCA; 2-60 keV) and the High Ener gy X-ray Timing Ex periment 



Table 1. RXTE PSR J0205 +6449 observation summary 



Obs. id. 


Date begin 


Date End 


MJD 


Exposure 17 
(ks) 


20259 


30-09-1997 


30-09-1997 


50721-50722 


16.96 


60130 


17-08-2001 


19-08-2001 


52138-52141 


80.35 


70089 


10-03-2002 


23-04-2003 


52343-52752 


268.95 


90080 


28-02-2004 


03-03-2005 


53063-53432 


243.23 


91063 


12-03-2005 


18-03-2006 


53441-53813 


328.74 



ting_t 

(HEXTE; 15-250 keV). The PCA dJahoda et all 1 19961) consists 



" Screened (GTI) exposure for PCA unit-2 



of five collimated Xenon proportional counter units (PCUs) with 
a total effective area of ~ 6500 cm 2 over a ~ 1° (FWHM) field of 
view. Each PCU has a front Propane anti-coincidence layer and 
three Xenon layers which provide the basic scientific data, and 
is sensitive to photons with energies in the range 2-60 keV. The 
energy resolution is about 18% at 6 keV. All data used in this 
work have been collected from observations in GoodXenon or 
GoodXenonwithPropane mode allowing high-time-resolution 
(0.9/js) studies in 256 spectral channels. 

The HEXTE instrument dRothschild et all 1 19981) consists of 
two independent detector clusters A and B, each containing four 
Na(Tl)/ CsI(Na) scintillation detectors. The HEXTE detectors 
are mechanically collimated to a ~ 1° (FWHM) field of view and 
cover the 15-250 keV energy range with an energy resolution of 
~ 15% at 60 keV. The collecting area is 1400 cm 2 taking into ac- 
count the loss of the spectral capabilities of one of the detectors. 
The best time resolution of the tagged events is 7.6/l/s. In its de- 
fault operation mode the field of view of each cluster is switched 
on and off source to provide instantaneous background measure- 
ments. Due to the co-alignment of HEXTE and the PCA, both 
instruments simultaneously observe the same field of view. 

RXTE observed PSR J0205 +6449 for the first time on Sept. 
30, 1997 (MJD 50 721) for abou t 17 ks . Data from this observa- 
tion were used by Murr ay et al.l d2002l) to confirm the pulsation 
discovered with Chandra. A dedicated much deeper observation 
was performed in the period August 17-19, 2001 (MJD 52138- 
52141), yielding about 80 ks good exposure time. Next, a moni- 
toring campaign started on March 10, 2003 (MJD 52343) which 
ended on April 23, 2003 (MJD 52752). The total (good) expo- 
sure time for this period was about 269 ks. A second monitoring 
round commenced on Feb. 28, 2004 and continued till March 
18, 2006 (MJD 53063-53813) yielding a total (good) exposure 
time of about 572 ks. A summary of all RXTE observations of 
PSR J0205 +6449 is given in Table Q] The total good-time expo- 
sure (after screening; see Sect. 13. Il l amounts 938.23 ks. 

2.2. XMM-Newton 

We searched the XMM-Newton observation database for obser- 
vations of the field around PSR J0 205 +6449 in which the EPIC- 
PN camera dStriider et all [2001) operated in Small- Window 
(SW) mode. This mode (4'A X 4'A field of view) offers suffi- 
cient time resolution (~ 5.67 ms) to sample the pulse-profile 
of PSRJ0205+6449 over the ~ 0.3-12 keV range. We found 
two observations (observation ids. 0004010101/0004010201) 
both performed on February 22, 2001 at 1'2 offset from 
PSRJ0205+6449 with durations of about 9.2 and 23.6 ks, re- 
spectively. 



L. Kuiper et al.: Hard X-ray characteristics of the energetic pulsar PSR J0205+6449 in 3C58 3 
Table 2. Phase-coherent ephemerides for PSR J0205 +6449 as derived from RXTE PC A (monitoring) data. 



Entry 


Start 


Purl 


to, Epoch 


V 


V 


V 


% 


Validity range 


# 


[MJD] 


[MJD] 


rMTD TDR1 


L±±z,j 


xlO" 11 Hz/s 


x10" 21 H7/S 2 









50721 


50722 


5U/2X.U 


15.230153381(92) 


-4.489 (nxea) 


U.U (nxea) 


0.2931 


2 


1 


52138 


52141 


52138.0 


15.224659060(16) 


-4.489(16) 


0.0 (fixed) 


0.5389 


4 


z. 






JZj'+J.V 






4-1 QftC\9k\ 






-x 

J 




jzj i J 










VI. 1 I J 1 


J 


4 1 


52639 


52752 


52639.0 


15.2227187742(4) 


-4.51063(1) 


0.0 (fixed) 


0.4147 


114 


5 


53063 


53173 


53063.0 


15.2211188258(7) 


-4.56381(2) 


0.0 (fixed) 


0.4814 


111 


6 


53173 


53312 


53173.0 


15.2206851395(13) 


-4.55932(4) 


+5.80(8) 


0.3170 


140 


7 


53312 


53401 


53312.0 


15.2201380229(35) 


-4.54739(19) 


+1.36(49) 


0.3458 


90 


8 


53401 


53469 


53401.0 


15.2197884405(37) 


-4.54444(25) 


+12.0(9) 


0.3273 


69 


9 


53469 


53546 


53469.0 


15.2195216777(26) 


-4.53041(15) 


+7.18(47) 


0.1970 


78 


10 


53546 


53637 


53546.0 


15.2192204400(23) 


-4.52131(12) 


+1.85(32) 


0.6030 


92 


11 


53637 


53726 


53637.0 


15.2188650817(24) 


-4.51951(12) 


+28.5(3) 


0.1413 


90 


12 


53726 


53813 


53726.0 


15.2185183927(26) 


-4.50300(14) 


+16.8(4) 


0.2436 


88 


13 2 


53726 


53764 


53745.0 


15.2184445054(20) 


-4.49904(15) 


-10.3(61) 


0.1147 


39 


14 2 


53750 


53814 


53782.0 


15.2183007045(13) 


-4.49502(8) 


+29.1(16) 


0.6055 


65 



A glitch occured between MJD 52515 and 52571 (see Livingstone et al., 2009, for more information). This entry describes 



the last part of the glitch recovery period. Its validity is questionable given the low number of TOA's, namely 4, translating 
to 1 degree of freedom in the TOA fit procedure. 

2 Ephemeris from Jodrell Bank radio data 

3 O is the phase offset to be applied to obtain consistent radio-alignment (see Equation[TJin Sect. l3.4t 



3. Timing 

3.1. RXTE PCA timing analysis 

The first step in the RXTE PCA data analysis was the screen- 
ing of the data. We generated good-time intervals (GTI) for each 
PCU separately, because the number of active PCU's at any in- 
stant was changing. Good time intervals have been determined 
for each PCU by including only time periods when the PCU in 
question is on, and during which the pointing direction is within 
0°05 from the target, the elevation angle above Earth's horizon 
is greater than 5°, a time delay of 30 minutes since the peak of a 
South-Atlantic-Anomaly passage holds, and a low background 
level due to contaminating electrons is observed. These good 
time intervals have subsequently been applied in the screening 
process to the data streams from each of the PCUs (e.g. see Table 
[T]for the resulting screened exposure of PCU-2 per observation 
run). 

Next, we selected event data from all three Xenon lay- 
ers of each PCU allowing us to better characterize the hard 
(> 10 keV) X-ray properties of PSR J0205+6449. The TT 
(Terrestial Time) arrival times of the selected events (for each 
sub-observation and for each PCU unit) have been converted 
to arrival times at the solar system barycenter (in TDB time 
scale) using 1) the JPL DE200 solar system ephemeris, 2) the in- 
stantaneous spacecraft position and 3) the sub-arcsecond celes- 
tial position of PSR J0205+6449. The position use d is: (a, 6) - 
(02 h 5 m 37 5 .92, +64°49'42'.'8) for epoch J2000 dSlane et all 
2002), which corresponds to (l,b)=(130.71931,3.08456) in 
Galactic coordinates. 



The pulse frequency and its first two time derivatives (v, v, v) 
were determined from PCA X-ray data solelfl demanding a 
maximum RMS value of only 0.01 period in the time-of-arrival 
(TOA) analysis. This requirement resulted in 13 timing models 
with validity intervals of typically 100 days. The ephemerides 
are listed in Table [2] In the TOA analysis we foll owed the steps 
outlined in Section 4 of iKuiper & Hermsenl (|2009l) . in this case, 
however, we made a high-statistics correlation template (show- 
ing clearly the two X-ray pulses) from the 80 ks observation 
during r un 60130. Our models ar e fully consistent with those de- 
rived bv lLivingstone et all d2009l) . who used a combination of X- 
ray (RXTE PCA) and radio (GBT and JBO) data. Also, we found 
evidence for the presence of two timing glitches using solely 
X-ray data, one occuring somewhere between MJD 52515 and 
52571 and a much stronger one occuring in t he RXTE monitor- 
ing g ap between MJD 52752 and 53063 (see lLivingstone et all 
2009, for more details on these glitches which they report to have 
fractional magnitudes Av/v ~ 3.4 x 10~ 7 and Av/v ~ 3.8 x 10" 6 , 
respectively). The frequency evolution history over the RXTE 
observation time stretch MJD 52138-53813 is shown in Fig.Q] 

The main difference between our work and that performed by 
[Livin gstone et al. (2009) is that we chose for an accurate (RMS 
< 0.01 ) description of the rotation behaviour of the pulsar with at 
most 3 timing parameters over a limited time stretch in stead of 
using many more timing parameters over a much wider time in- 
terval. In the latter approach the need for the (unphysical) higher 
order timing parameters reflects the presence of (strong) timing 
noise. 



3.2. Timing solutions: ephemerides 

We generated pulsar timing models (ephemerides) specifying 
the rotation behaviour of the pulsar over a certain time stretch. 



1 X-ray timing data are not hampered by time-variable dispersion 
measure (DM) variations as is the case for the radio data, and there- 
fore the remaining scatter in the obtained timing solutions is less. This 
has been verified for the Crab pulsar. 
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Fig. 1. Evolution of the spin-frequency of PSR J0205 +6449 as 
derived from RXTE PCA data over the period MJD 52138- 
53813 with respect to the linear trend of the phase coherent tim- 
ing model of period MJD 52433-52515 (entry #3 of Table 01. 
As solid lines the entries #1-12 of Table [2] are plotted, while 
for the period 52571-52752 also the linear fit to the incoher- 
ent frequency measurements is shown as solid line. Incoherent 
frequency measurements over the period MJD 52544-52752 are 
shown as data points. Note the presence of (at least) two timing 
glitches: one in the period 52515 to 52571 and a second stronger 
one between 52752 and 53063 (see text). 



3.3. X-ray/radio pulse profile phase alignment 

The Jodrell Bank observatory (JBO) made observations at a ra- 
dio frequency of 1.4 GHz from MJD 53725 to 54666, and there- 
fore overlaps for about 89 days with the second RXTE monitor- 
ing cycle in the period MJD 53725 to 53813. For two time seg- 
ments in this interval, MJD 53726-53764 (number of TOAs, 28) 
and MJD 53750-53814 (number of TOAs, 29), accurate (RMS 
< 0.01) timing models are constructed with 3 timing parame- 
ters (see also Table |2J. The 1.4 GHz single-pulse radio profile 
(in 400 bins) is shown in Fig. [3^ with a fiducial point (defin- 
ing radio-phase 0.0) corresponding to the centre of gravity of 
the single pulse (just before the pulse maximum). Folding the 
barycentered X-ray time tags from period MJD 53726 to 53813 
upon these radio-ephemerides put the main X-ray pulse (pulse- 
1) at phase 0.089 + 0.001 (statistical error onl y), consistent with 
the va lue quoted for the JBO-PCA offset in Livingstone et al.l 
(2009), namely 0.085 + 0.010. Next, we determine through cor- 
relation analysis the phase shifts to be applied to the X-ray pulse 
profiles from the data periods of entries 0-12 of Table|2]to align 
these to the radio-aligned X-ray profile of period MJD 53726- 
53813. These shifts (<t>o) are given in Table|2] 



1.85x10' 




Fig. 2. High-statistics radio-aligned PCA pulse phase distribu- 
tion in 180 bins for PHA range 5-44 (~2-20 keV) combining 
all available radio-aligned pulse phase distributions from differ- 
ent data segments. Error bars represent lcr uncertainties. The 
best fit (x 2 = 199.47 for 171 degrees of freedom) model com- 
posed of two asymmetric Lorentzians plus background is super- 
posed as dashed red line. The rising wings of the pulses are 
steeper than the trailing wings. 

as shown in Table [2] Proper X-ray/radio phase alignment, <£>(?), 
is obtained by subtracting <t>o as shown in the following formula: 



O(f) 



•ft- 



1 



1 



t ) + -v-(t-t y + -v-(t- 

l 6 



to) - $0 



(1) 



Combining the radio-aligned phase information for all PCA- 
data covered with a proper ephemeris (see Table |2]l yielded an 
event matrix N(0,E) of (180x256) elements. For this purpose 
we binned the pulse-phase interval [0,1] into 180 phase bins for 
all 256 PCA PHA channels. From this matrix a high-statistics 
pulse-phase distribution was extracted for the PHA channels 5- 
44 (~ 2-20 keV). This distribution is shown in Fig. [2] 

3.5. X-ray pulse profile characterization 

Initially, we fitted, analogous to Livingstone et al. (2009), our 
high-statistics RXTE PCA pulse profile shown in Fig. [2] with 
a model consisting of 2 Gaussians, each with free scale, width 
and position, plus background. However, this model rendered a 
poor/unacceptable fit (x 1 = 265.48 for 180-7 degrees of free- 
dom). Next, we tried a double symmetric Lorentzian model plus 
background in order to give more weight to the wings of the 
pulses. This model provided a better description of the measured 
pulse-phase distribution (y 2 = 239.13 for 180-7 degrees of 
freedom), but is still poor. Finally, we abandoned the descrip- 
tion in terms of symmetric functions and used a combination of 
2 asymmetric Lorentzians plus background. This model (9 free 
parameters) is specified below: 



3.4. Combined X-ray event matrix from PCA observations B,pi,p 2 ) = B + N\(<p;p{) + Niifcpz) 



(2) 



Barycentered PCA X-ray event times are finally folded upon an 
appropriate timing model composed of v, v, v and the epoch fo, 



In this formula p\ represents the 4 model parameters, 
(N\,$\,Y\uY\r), of the first asymmetric Lorentzian, pz the 
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Table 3. X-ray pulse profile characterization of 
PSR J0205 +6449 from a fit involving two asymmetric 
Lorentzians plus background 



Parameter 



Value 



lir-error 



Pulse- 1 





0.0831 ± 


0.0004 


r„ 


0.0175 ± 


0.0009 




0.0252 ± 


0.0011 




Pulse-2 






0.5709 ± 


0.0015 


r 2; 


0.0214 ± 


0.0036 


r 2r 


0.0502 ± 


0.0056 




Derived quantities 


<D 2 - 0[ 


0.488 ± 


0.002 


Ni/N 2 


3.72 ± 


0.23 


n 


0.688 ± 


0.011 


h 


0.312 ± 


0.014 


R = hlh 


2.2 ± 


0.1 


n = (t u 


+ T lr )/2 0.0214 ± 


0.0007 




1.41 ± 


0.05 ms 




+ r 2 ,.)/2 0.0358 ± 


0.0034 




2.35 ± 


0.22 ms 



" Statistical error only, the systematic error is of the 
order of 0.01 (see Livingstone et al., 2009) 

* Relative contribution of the integrated flux in pulse- 1 
to the total pulsed flux 



JBO: Radio 





RXTE PCA: X-rays 




FERMI LAT: y-rays 



0.0 0.5 



1.0 



1.5 



2.0 



equivalent parameters, (N2, <Pi, ^21, ^2r), describing the second 
asymmetric Lorentzian and B is the value of the background 
level. The first asymmetric Lorentzian is described by the fol- 
lowing expression: 



Ni(4>;pi) 



Ni 



U-<Pi)l(Tul2)Y + \ 

((<M>,)/(r lr /2)) 2 +i 



(3) 



In this description N\ is the maximum value of pulse- 1 
reached at <j>\, the location of the maximum of pulse-1, Tu/2 
is the width of the left wing of the pulse-1, and finally Ti r /2 
is the width of the right wing of the pulse-1. A similar expres- 
sion and equivalent definitions hold for the second asymmetric 
Lorentzian. 

This composite model provided an excellent fit, x 1 = 199.47 
for 171 degrees of freedom, with best-fit parameters and their 
lcr error estimates listed in Table [3] (see also the best-fit model 
superposed on the data in Fig- EJ. A description in terms of two 
asymmetric Lorentzians plus background provides a 7.7cr im- 
provement over the two-Gaussians-plus-background model and 
a 6cr improvement over the two-Lorentzians-plus-background 
model, taking into account the two (=9-7) additional degrees of 
freedom in both cases. Therefore, our analysis do es not support 
the assumption made bv lLivingstone et all (12009) of an underly- 
ing double Gaussian shape for the X-ray profile. We find the X- 
ray pulses to be sharper, especially for pulse-2. For both pulses 
the rising wings are significantly steeper than the trailing wings. 

The X-ray peak separation <p2 - <p\ derived in this work 
is 0.488(2), significantly smaller than the value estimated by 



Pulse Phase 

Fig. 3. A comparison in absolute phase of the JBO radio (1.4 
GHz), RXTE-PCA (-2-20 keV) and Fermi LAT (> 100 MeV) 
pulse profiles. Note the change in the relative contributions of PI 
and P2 in the X-ray and y-ray windows. 



Livingsto ne et alJ d2009l) . but consistent with the s eparation of 
0.49+0 .01 ±0.0lB measured at high-energy y-rays bv lAbdo et al.l 
(2009a) using Fermi LAT > 100 MeV data. The comparison of 
the shapes and absolute phases of the JBO radio, our RXTE- 
PCA X-ray and the Fermi-LAT profiles is shown in Fig. [3] The 
main X-ray pulse (PI) appears to be the sharpest pulse in this 
comparison. 



3.6. X-ray pulse profile variability 

We investigated the stability of the X-ray pulse-shape as a func- 
tion of time. Therefore, we fitted the measured X-ray pulse- 
phase distribution (PHA range [4,27] ~ 2-11 keV) for 15 time 
periods in terms of a constant background and the shapes of 
pulse-1 and pulse-2, separately. The times of these data points 
correspond to those of the X-ray timing models shown in Table 
[2] (entries to 11; 12 points), augmented with two measure- 
ments during the last RXTE monitoring period, MJD 53726- 
53813, covering entry #12 and finally, with a data point cover- 
ing the post-glitch- 1 period MJD 52544-52607, yielding eventu- 
ally 15 independent measurements. The splitting of period MJD 



2 The first error s pecifies the statisti cal error and the second the sys- 
tematical error (see lAbdo etaLll2009aL for more details) 
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Fig. 4. The ratio of the integrated flux in pulse- 1 over pulse-2 
as a function of time for the PCA energy band ~ 2-1 1 keV. One 
data point at time interval MJD 53736-53749 (i.e. the "anomaly" 
period) deviates ~ 3.5cr (single trial) from the time-averaged 
value of 2.2 ±0.1. 



53726-538 13 into the intervals MJD 53736-53749 (2 RXTE sub- 
observations) and MJD 53760-53813 (5 sub-observations) was 
driven by the detection of a "timing anomaly" in the former in- 
terval during the phase-coherent timing analysis. At a later stage 
of this work it turned out that this "anomaly" was caused by in- 
correct RXTE clock corrections just after the introduction of a 
leap second on 2006, January 1 . 

The profile fitting procedure yields the flux ratio R = 
(see for the definition Tabled for each time interval. The results, 
R(t) vs. t, are shown in Fig. [4] with superposed as long-dashed 
line the Pl/P2-flux ratio from the time-averaged high-statistics 
profile (R = 2.2±0.1; see Tabled along with its 1 -<x error region 
(shaded area). One data point, corresponding to the "anomaly" 
period, deviates ~ 3.5<x from the time-averaged value. Taking 
into account the number of trials (15) its significance reduces 
to 2.7<x, still indicating an interesting hint for variability. The 
pulse-phase distribution during the "anomaly" period is shown 
in Fig. [5] In this figure we also superposed as dotted line the best 
fit model in which the shapes (two asymmetric Lorentzians) for 
each of the two pulses are identical to those derived for the time- 
averaged profile in Fig. [2] and detailed in Table [3] 

Compared to the other measurements, where P2 is some- 
times hardly visible, we see strongly enhanced P2 emission dur- 
ing this period. From spectral analysis of the PI and P2 emis- 
sions during the "anomaly" period it turns out that the PI flux 
is comparable to its time-averaged value contrary to the P2 flux, 
which shows a clear enhancement by almost a factor 2. 

This leads to the conclusion that we see an interesting indi- 
cation for flux variability for P2 without a change of its pulse 
shape. Finally, we checked the JBO radio profile assembled dur- 
ing the "anomaly" period for possible morphology changes e.g. 
the appearance of a new feature, but we found none. 

3. 7. RXTE HEXTE timing analysis 

HEXTE operated in its default rocking mode during the observa- 
tions listed in Table [T] allowing the collection of real-time back- 
ground data from two independent positions ± 1°5 to either side 
of the on-source position. For the timing analysis we selected 
only the on-source data from both clusters. Good-time inter- 
vals have been determined using similar screening filters as used 
in the case of the PCA. The selected on-source HEXTE event 
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Fig. 5. The pulse profile of PSR J0205+6449 (60 bins) in the ~2- 
1 1 keV band during the "anomaly" interval (MJD 53736-53749; 
2 RXTE sub-observations). Strongly enhanced P2 emission is 
detected. The best fit model, composed of a background plus two 
asymmetric Lorentzians of the same shapes as shown in Fig. 2, 
is superposed as dashed line. 



times have subsequently been barycentered and folded upon the 
ephemerides listed in Table [2] taking into account proper radio- 
phase referencing. Thus, we obtained time-averaged HEXTE 
pulse phase distributions in 256 spectral channels (15 - 250 keV) 
for the combination of observations listed in Table Q] The total 
dead-time corrected exposure time collected for clusters A and 
B amounts, 400.6 ks and 426.3 ks, respectively. Pulse profiles 
for the bandfl 14.7-28 and 33. 1-132.6 keV, are shown in panels 
c and d of Fig. [6] Fitting a model, comprising the (asymmetric) 
Lorentzians shapes of Pulses 1 and 2 and a flat background, to 
these phase distributions yielded detection significances of 7.8<x 
and 3.8<xfor the 14.7-28 and 33. 1-132.6 keV bands, respectively 
(2.9cr for the band, 64. 1-132.6 keV). Therefore, pulsed emission 
of PSR J0205 +6449 has been detected up to ~ 132 keV, well 
above the sensitivity band of the PCA. 

3.8. XMM-Newton timing analysis 

The XMM EPIC-PN data were screened for solar (soft proton) 
flares by creating a light curve for events with energies in ex- 
cess of 10 keV. From the resulting count rate distribution, as- 
sumed to be Gaussian in absence of any flares, we could iden- 
tify periods during which the rate exceeds its mean value plus 
three times the width of the distribution. These periods are ig- 
nored in subsequent analysis. Next, we selected events from a 
sufficiently large circular region centered on PSR J0205+6449 
with a radius of 60" to ensure that all pulsar counts are included 
and barycentered the event times of these events. Because the 
XMM-Newton observations have been performed before the 80 
ks RXTE observation (60130) no valid ephemeris was available 
for the XMM data period. Therefore, we performed a limited pe- 
riodicity search around the predicted frequency value based on 



3 We ignored spectral data from the band 28-33. 1 keV because of the 
presence of a huge background line originating from the activation of 
Iodine. 
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ous to our X-ray pulse profile variability study shown in Sect. 
13.61 ) to the measured pulse-phase distribution, ^(O), in various 
user-selected energy bands: 
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Fig. 6. XMM-Newton EPIC-PN pulse profiles (30 bins) of 
PSRJ0205+6449 for the 0.5-3 and 3-12 keV energy bands 
(panels a,b). Panels c and d show the RXTE HEXTE profiles 
(60 bins) for the 14.7-28 and 33.1-132.8 keV energy bands. 
Significant pulsed emission is detected up to ~ 130 keV and 
down to ~ 0.95 keV. 

entry- 1 of Table [2] We found a ~ 10cr signal right at the pre- 
dicted frequency using only events with energies between 4 and 
12 keV. The folded pulse profile was compatible with the high- 
statistics PCA profile (see Fig.|2]i taking into account the blurring 
of the pulse profile due to the limited time resolution of 5.67 ms 
(=0.086 in phase units). For subsequent studies we used an ex- 
traction radius of 15" because the signal-to-noise ratio peaks at 
that value. We created pulse-phase distributions, combining both 
adjacent XMM observations, in 30 bins (i.e. oversampled by a 
factor of ~ 3) for 585 energy intervals over the 0.3 to 12 keV 
range, each 0.02 keV wide. The (radio-aligned) EPIC-PN pulse 
profiles for the 0.5-3 and 3-12 keV bands are shown in the panels 
a and b of Fig. [6] an d have signifi cances of 5.8 and 9.7<x (adopt- 
ing Z 2 -test; Bucche ri et allfl9 83). respectively. Pulsed emission 
has been detected down to ~ 0.95 keV. 



4. Pulsed spectra from RXTE PCA/HEXTE and 
XMM-Newton data: Total pulsed, P1 and P2 

From the pulse-phase distributions N(<$>, E), derived for RXTE 
PCA, HEXTE and XMM-Newton EPIC PN, we extracted pulsed 
excess counts by fitting the following model function (analogu- 



Af(<D) = b + c 1 xT i (<Z>) + c 2 x r 2 (0) 



(4) 



In this formula b represents the (constant) unpulsed/DC level, c\ 
and C2 the scales of the two asymmetric Lorentzian templates, T\ 
and T 2 (both normalized to 1), respectively (see Sect. 13.5b . This 
model provided statistically good fits to all PCA and HEXTE 
profiles. Good fits could be derived for the EPIC-PN profiles 
after convolving the Lorentzian templates with the poorer time 
resolution. 

For each instrument the pulsed excess counts in the various 
energy bands for the first (PI) and the second pulse (P2) and 
the sum (=total pulsed, TP) can be translated to photon fluxes 
provided that proper energy response matrices are used. 

In the case of the PCA we constructed time-averaged en- 
ergy response matrices for each PCU separately taking into ac- 
count the different (screened) exposure times of the involved 
PCU's during the time period of interest. For this purpose we 
used the ftools version 6.4 programs pcarsp and addrmf. To con- 
vert PHA channels to measured energy values, £pha, for PCU 
combined/stacked products we also generated a weighted PCU- 
combined energy response matrix. 

For HEXTE we employed cluster A and B energy-response 
matrices separately, taking into account the different screened 
on-source exposure times and the reduction in efficiency in case 
of off-axis observations. The on-source exposure times for both 
clusters have been corrected for considerable dead-time effects. 

Finally, we created energy response files (effective area {arf) 
and energy redistribution matrix (rmf)) for the EPIC PN operat- 
ing in small window mode taking into account the reduction in 
effective area given the 15" source extraction radius used. For 
this purpose we employed the XMM SAS (vrs. 7.1.0) software 
tools arfgen 1.73.3 and rmfgen 1.55.1. 

We assume simple power-law models in the form, F y 
K ■ {E y IEoY r with F the photon-index and K the normaliza- 
tion in ph/cm 2 s keV at the pivot energy E$, for the underlying 
photon spectra of PI, P2 and its sum TP. We fixed the absorb- 
ing interstellar Hydrogen column Nh to 3.4 x 10 2 1 cm" 2 (see 
the "P L-model for neutron star" entry in Table 2 of ISlane et all 
2004). These models have been fitted in a forward folding pro- 
cedure using the appropriate response matrices to obtain the op- 
timum spectral parameters, K and F, and the reconstructed spec- 
tral flux points from the observed pulsed count rates. We veri- 
fied that the measured high-statistics RXTE-PCA spectrum, as 
well as the EPIC-PN spectrum and the total spectrum includ- 
ing also HEXTE data are fully consistent with this non-thermal 
simple power-law model. There is in the pulsed X-ray spectrum 
above ~ 0.5 keV no indication for a thermal black-body compo- 
nent, a conclusion also reached for the tota l emission f r om th e 
compact source by ISlane et al] d2002l) and ISlane et al] (|2004), 
who reported a powe r -law s pectral index of ~ 1.7. We note that 
Kargaltse v & Pavlovl ((2008) in their Table I erroneously mark 
this pulsar to have a black-body component. Furthermore, in this 
work we only show the unabsorbed spectra i.e. the interstellar 
absorption has been modeled out. 

In Table|4]the best fit values are listed of the spectral param- 
eters for the total pulsed emission TP, and emissions of PI and 
P2 using PCA data only, and for TP using the EPIC PN, PCA 
and HEXTE combination over the extended energy band 0.56 to 
267.5 keV. All spectra have a consistent shape with index ~ 1 .03. 
We note for comparison, that for energies > 100 MeV the Fermi 
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Table 4. Best fit values for the photon flux spectra of the total 
pulsed emission (TP), the first (PI) and second (P2) pulse emis- 
sions of PSR J0205 +6449 assuming a power-law model of the 
form F y = K ■ (E 7 /E Q )- r . 



Parameter 



TP 



PI 



P2 



PCA (2.5-54.0 keV) 

#(l(T 6 ph/cm 2 skeV) 2.93 ± 0.05 2.03 ± 0.03 0.90 ± 0.04 

T 1.06 ±0.03 1.04 ±0.03 1.10 ±0.08 

E (keV) 8.34 8.34 8.34 

F2-30 (10- 13 erg/cm 2 s) 10.67 ±0.16 7.43 ±0.11 3.24 ±0.13 

XMM/PCA/HEXTE (0.56-267.5 keV) 

K (10~ 6 ph/cm 2 skeV) 2.85 ± 0.04 

T 1.03 ± 0.02 

E (keV) 8.49 

Fo.5-i50 (10~ 12 erg/cm 2 s) 5.48 ± 0.28 



LAT measures for TP a much softe r spectrum w ith spectral index 
2.1 with cutoff at ~ 3 GeV dAbdo et all l2009al) . and that P2 ex- 
hibits at high-energy y-rays a significantly harder spectrum than 
PI. 

The photon spectrum (vF v representation) over the 0.56 - 
267.5 keV energy band of the total pulsed emission combin- 
ing XMM-Newton EPIC-PN, RXTE-PCA and HEXTE data, as 
derived in this work, is shown in Fig. Q in a much wider en- 
ergy frame (0.1 keV-10 GeV) by including the best fit and its 
uncertainty r ange m easured by Fermi for energies > 100 MeV 
dAbdo et al.L |2009a). The luminosity of the pulsed emission of 
PSR J0205 +6449 apparently reaches a maximum in the MeV 
band. For comparison are also shown the total pulsed emission 
spectra of the Crab, PSR B1509-58 as well as the "middle-aged" 
Vela pulsar. 



5. Summary 

In this paper we derived for the young rotation-powered pulsar 
PSR J0205 +6449 the timing and spectral characteristics over the 
broad X-ray band from ~ 0.5 to ~ 270 keV, using data from the 
RXTE PCA and HEXTE, and XMM-Newton EPIC PN. These 
X-ray characteristics complement our knowledge about this pul- 
sar in the radio domain and the high-energy y-ray band for ener- 
gies above 100 MeV. 

- Our phase-coherent ep hemerides (see Table [2]) are con- 
sistent with those derived by Livingst one et al.l d2009l) with the 
main difference that we used solely X-ray data (from the RXTE 
PCA) and fitted at most three timing parameters (v, v, v) over 
more limited time intervals. 

- The X-ray pulse profile consists of two sharp pulses which 
can be described with 2 asymmetric Lorentzians, each with the 
rising wing steeper than the trailing wing, and full-width-half- 
maximum 1.41 ± 0.05 ms and 2.35 + 0.22 ms, respectively. 
These profiles are sharper than reported by [Livingstone et al.l 
(120091) . 

- The first X-ray pulse lags the single radio pulse in phase 
by 0.089 + 0.001 (statistical error); the phase separation be- 
tween the two X-ray pulses amounts 0.488 + 0.002, fully con- 
sistent with the value 0.49 ± 0.01 ± 0.01 (statistical and sys- 



tematic errors) repo rted for high-energy y-rays above 100 MeV 
dAbdo et aUl2009ah . 

- We find an indication for a flux increase by a factor ~ 2, 
~ 3.5cr above the time-averaged value, for the second, weaker 
pulse during a two-week time interval, while its pulse shape did 
not change. During this time window, the morphology of the 
JBO radio profile of PSR J0205 +6449 did not change, notably, 
there was no indication for a second pulse. 

- We detected the pulsed signal significantly for the first time 
down to ~ 0.95 keV with XMM-Newton EPIC PN, and up to ~ 
130 keV by analysing RXTE HEXTE data. The morphologies of 
the EPIC PN (taking into account the coarser timing resolution) 
and the HEXTE profiles are consistent with that measured with 
the PCA. 

- The spectrum of the pulsed X-ray emission is of non- 
thermal origin, exhibiting a power-law shape with photon index 
r = 1.06 ± 0.03, fitting just the high-statistics PCA data, and T 
= 1.03 ± 0.02, fitting over the broader energy band from ~ 0.5 
to ~ 270 keV by including also the EPIC-PN and HEXTE flux 
values. There is no indication for a black-body component in the 
soft X-ray spectrum above 0.5 keV. 

- We do not see a spectral difference between the spectra of 
the two X-ray pulses in the PCA data. Both spectral photon in- 
dices are fully consistent with the time averaged value for the to- 
tal pulsed emission (see Table|4]i. Note that the relative strengths 
of PI and P2 in the X-ray and high-energy y-ray windows re- 
verse (see Fig. [3}; the spectrum of P2 has to harden significantly 
with respect to that of PI between a few hundred keV and 100 
MeV 



6. Discussion and Conclusions 

In the introduction we noted that PSR J0205 +6449 is now one 
of only three young (< 10,000 year old) pulsars which are de- 
tected in the classical X-ray band and at hard X-rays above 20 
keV, as well as at high-energy (> 100 MeV) y-rays, the oth- 
ers being the Crab pulsar and PSR B 1509-58. Fig. [7|shows that 
these three young pulsars reach their maximum luminosities be- 
low 100 MeV, while the "middle-aged" Vela pulsar (character- 
istic age 1 1 .4 kyr) reaches its maximum at GeV energies (for 
the latest results on the Vela pulsar for ener gies above 100 MeV, 
see the Fermi results by Abd o et al.Ll200 9b). The latter spectrum 
is characteristic for older pulsars reported to be detected above 
100 MeV (e.g. see t he first Fermi Larg e Area Telescope catalog 
of y-ray pulsars by I Abdo et al.l l2009d) 

Comparing in more detail the high-energy spectra of the 
young pulsars in Figure [7] we notice large differences. For en- 
ergies below 10 keV the flux values of PSR J0205 +6449 are ~ 
4 orders of magnitude below those of the Crab, while around 
10 MeV the difference is reduced to about a factor of 10. 
The X-ray spectrum of PSR J0205 +6449 is, thus, very much 
harder than that of the Crab. The total high-energy spectrum 
of PSR J0205 +6449 appears to reach its maximum luminos- 
ity at MeV energies, like is the case for PSR B1509-58/PSR 
J1513-5908. The spectral break for the latter spectrum between 
10 MeV and 100 MeV (see flux values in Figure [7]l measured 
with COMPTEL a nd EGRET aboard the Compton Gamma- 
Ray Observatory bv | Kuiper et al.l d 19991) has been confirmed by 
iPellizzoni etafl d2009l) . These authors report for PSR B1509-58 
a softening of the photon index T from ~ 1.7 to ~2.5 going from 
tens to hundreds of MeV (but do not provide pulsed-flux values). 

Interestingly, the X-ray spectrum above 2 keV of 
PSR J0205 +6449 resembles that of the slightly older Vela pul- 
sar (similar spectral index), but the Lx/L y ratio for the pulsed 
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Table 5. Characteristics of the three young (< 10 kyr) X-ray and y-ray emitting pulsars in comparison with the middle-aged Vela 
pulsar (PSR B0833-045). The luminosities of the pulsed emission L are calculated as L — AncpFfa with values for the distance d 
taken from the table, and the beaming fraction /n set to 1 . F represent the pulsed flux. 



Source 


d 


P 


r 


L sc i 










K 






4 


4 




[kpc] 


[ms] 


[ky] 


[erg/s] 


[erg/cm 2 s] 


[erj 


|/cm s] 




[erg/s] 






[erg/s] 




B0531+21 


2.0 


29.7 


1.3 


4.4E+38 


(5.68 ± 0.05)E-09 


(1.3 ± 


0.1)E-09 


(2.72 


± 0.02)E+36 


6.2E-3 


(6.1 


± 0.3)E+35 


1.4E-3 


B 1509-58 


5.8 


151.5 


1.6 


1.7E+37 


(1.46±0.02)E-10 


(5.1 ± 


2.5)E-11 C 


(5.88 


± 0.08)E+35 


3.5E-2 


(2.1 


± 1.0)E+35 r 


1.2E-2 C 


J0205+6449 


3.2 


65.7 


5.4 


2.7E+37 


(0.36±0.02)E-11 


(6.7 + 


0.5)E-11 


(4.45 


± 0.20)E+33 


1.7E-4 


(8.2 


± 0.6)E+34 


3.0E-3 


B0833-045 


0.287 


89.3 


11.4 


6.9E+36 


(0.88 + 0.29)E-ll 


(7.9 + 


0.3)E-09 


(8.67 


±2.85)E+31 


1.3E-5 


(7.8 


± 0.3)E+34 


1.1E-2 



" Luminosities, fluxes and efficiencies labeled with x are evaluated for the 2-100 keV band. 
h Luminosities, fluxes and efficiencies labeled with y are evaluated for the 0.1-10 GeV band. 

c The y-ray energy flux of PSR B1509-58 in the 0.1-10 GeV band has been derived from the (total) photon flux values for the 100-300 and 
300-1000 MeV bands as given in Kuiper et al. (1999) assuming a power-law shape with photon index of 2.5, and should be considered as 
an upper-limit to the pulsed flux of PSR B1509-58 in the 0.1-10 GeV band. 



component differs by a factor -50; higher for PSR J0205 +6449 
(see Table [5] which is introduced below). The LxlL y ratio of 
PSR J0205 +6449 is in between those of Vela and PSR B1509- 
58, namely, the L x /L y ratio for PSR J0205 +6449 is a factor 
-50 smaller than that for PSR B 1509-58. Note, that for the 
quoted flux ratios it is assumed that the beaming fractions in 
the X-ray and y-ray bands are the same. We know, however, that 
these are in many cases different. More importantly, in the X- 
ray spectra below e.g. 2 keV there are no indications for black- 
body components in the spectra of Crab, PSR B 1509-58 and 
PSRJ0205+6449. In contrast, the (pulsed) Vela spectrum ex- 
hibit s below 2 keV a bl ack-body peak (not shown in Fig. [71 see 
e.g. lPavlov etaTl 120011) . which is characteristic for middle-aged 
and older rotation powered pulsars. Therefore, the spectral prop- 
erties of PSR J0205 +6449 confirm that we are dealing with a 
young pulsar, and suggest a real age between those of Vela and 
PSR B 1509-58, favouring its characteristic age of 5.4 kyr over 
that of SN 1181 (828 yr). 

Table [5] lists for the four pulsars discussed above in order of 
characteristic age (r = P/2P) the spin-down luminosities 
and fluxes F and luminosities L in the X-ray 2-100 keV and 
gamma-ray 0. 1-10 GeV bands, as well as the corresponding effi- 
ciencies to convert spin-down energy into emission in these en- 
ergy bands. The luminosities are calculated as L = 4nd 2 Ffn, 
with the values for the distance d taken from the table, and the 
value for /n, which is the beaming fraction, set to 1 . At first sight 
one could argue that there is an evident anticorrelation between 
characteristic age and X-ray luminosity, independent from dif- 
ferences in the beaming fractions, but this becomes less obvious 
when we consider the X-ray efficiencies instead of luminosities. 
In the gamma-ray band there is not any indication for a correla- 
tion. The listed gamma-ray efficiencies differ less than a factor 
~10, ignoring differences in beaming fraction, while the latter 
differences can be substantial. 

There are also large differences in the morphologies of the 
pulse profiles of the three young pulsars. Comparing the pulse 
profiles detected for Crab and PSR J0205 +6449 at X-ray ener- 
gies and high-energy y-rays, then there are also some simulari- 
ties: both exhibit two pulses with peaks separated ~0.5 and ~0.4 
in pulse phase, respectively, and the X-ray and y-ray pulses are 
aligned in phase. However, the pulses in the Crab profile are sig- 
nificantly broader than those of PSR J0205 +6449 and emission 
is also detected between the two Crab pulses. The latter is not the 
case for the X-ray profile of PSR J0205+6449, but, interestingly, 




10~ 4 10" 2 10° 10 2 10* 



E [MeV] 

Fig. 7. Broad-band total-pulsed photon spectrum of 
PSR J0205+6449 (aqua colored) compared with the pulsed 
spectra of PSR B0531+21 (Crab; red), PSR B0833-45 (Vela; 
dark blue) and PSR B 1509-58 (purple). The (hard) X-ray 
spectrum (0.56-267.5 keV) of PSR J0205 +6449 has been 
derived in this work, and the best-fit power-law model (index 
~ 1.03) has been superposed. The > 100 MeV spectrum of 
PSR J0205 +6449 is the model fit to the Fermi spectrum from 
lAbdo et all ([2009a). Also for Crab and Vela the best-model 
fits to the recently p ublished Fermi spec tra for energies > 100 
MeV are shown (see lAbdo et all l2009blldl respectively). For all 
spectra the interstellar absorption has been modeled out (only 
effective below ~ 5 keV). 



emission between the pulses has been detected in the Fermi pro- 
file at the 5cr level. 

Furthermore, the Crab main radio pulse is in general phase 
coincidence with the broad X-ray/y-ray pulse. The peak of this 
main radio pulse lags that of the first X-ray/y-ray pulse in phase 
by only ~ 0.008 or 280 /zs; see for consistent estim ates from 
INTE GRAL, RXT E and EGRET iKuiperetaD (jool, and from 
Fermi lAbdo et al.l (l2009dl) . On the other hand, the Crab radio 
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precursor preceeds the first, main X-ray/y-ray pulse in phase by 
~0.09, or 3.2 ms, being located around the start of the leading 
wing of the high-energy pulse. In the case of PSR J0205 +6449, 
the single radio pulse is also preceeding the first narrow X-ray/y- 
ray pulse in phase by ~ 0.083 or 5.4 ms, and is fully separated 
in phase, the radio pulse being located just before the onset 
of the first high-energy pulse (see Fig. 0). This strongly sug- 
gests that the analogue of the radio pulse of PSR J0205 +6449 
is the weak radio precursor of the Crab. Also, that there are 
no counterparts in the radio profile of PSR J0205 +6449 to the 
two high-energy pulses of PSR J0205 +6449, contrary to the sit- 
uation for the Crab. This means that for this young pulsar ex- 
hibiting sharp non-thermal high-energy pulses, we do not see 
evidence for radio emission originating from the same site in 
the magnetosphere, e .g. in slot gaps (two-pole caustic emission, 
iDvks & RudakLl2003h or ou ter gaps ( outer- magnetosphere emis- 
sion. see lCheng et al.L[T9l6llRomanil[T996l;lHirotanill2006l from 
a region close to the light cylinder). Possibly, this radio com- 
ponent of PSR J0205 +6449 is just too weak to be detectable, 
but might be revealed in a search for giant radio pulses in the 
phase intervals of the high-energy pulses. Namely, for a number 
of young and milli-second radio pulsars phase coincidences be- 
tween the high-energy pulses and giant radio pulses have been 
reported. Two examples: the Crab for which the distribution of 
giant radio pulses is remarkably similar to the average em ission 
profile of the radio main and interpulse dPopov e t all 2006) and 
milli- second pulsar PSR B1937+21 for which [Cusumano et all 
(2003) reported the phase coincidence of two sharp high-energy 
X-ray pulses with two phase intervals exhibiting giant radio 
pulses, which trail the two normal radio pulses. The latter ex- 
ample might be revealed for PSR J0205+6449. 

The high-energy pulse profile of PSR B1509-58 differs to- 
tally from those of Crab and PSR J0205+6449. At hard X-rays 
and soft y-ray s below 10 MeV the profile consists of a single 
structured broad pulse, which can be explained as being com- 
posed of two Gaussian pulse profiles s eparat e d ~0.14 in phase 
with different spectra (Ku iper et all 1 1999b ICusumano et all 
1200 lb . the second broader pulse peaking at ~0.35, with the main 
radio pulse at phase 0. Above 10 MeV, the COMPTEL pro- 
file between 10 and 30 MeV and the EGRET profile between 
30 and 100 MeV suggest the presence of a n additional high- 
energy pulse at phase ~0.85 (Kuipe r et all I1999I) . The latter 
se ems now to be confirme d in the AGILE profile of PSR B1509- 
58 lPellizzoni et al. ( 2009), which shows the main pulse for ener- 
gies above 100 MeV at phase ~0.35, and a second possible pulse 
at ~0.85. It is now ambiguous what phase difference between 
high-energy pulses of PSR B1509-58 (~0.14or -0.5) should be 
considered for comparison with the morphology of pulse profiles 
of the other young pulsars. 

The above cited different models aiming to explain the pro- 
duction of non-thermal high-energy emission in the magneto- 
spheres of rotation-powered pulsars do not address flux vari- 
ability. Furthermore, there was also no observational evidence 
for such variability till the ma gnetar-like outburst of the high- 
field pulsar PSR J1846-0258 (iGavriil et all 120081) . which de- 
cayed with an 1/e-t i me c onstant of ~ 55 days. It was shown by 
Kuip er & Hermsenl (120091) that the radiative outburst was trig- 
gered by a major spin-up glitch, and that, most interestingly, 
the shape of the X-ray pulse profile did not change during the 
outburst. For the flux increase by a factor of ~ 2 of the non- 
thermal emission from the second pulse of PSR J0205 +6449 
during a two-week time period, we did not see a variation in 
pulse shape, either. However, there was no indication for glitch- 
ing activity. The significance of the variability is insufficient to 



draw strong conclusions, but it seems waranted to start searching 
for such variability in the emission from the increasing sample 
of rotation-powered pulsars emitting non-thermal high-energy 
emission. 

In conclusion, we accurately measured for the young rotation 
powered pulsar PSR J0205+6449 the morphology of the X-ray 
light curve and the spectrum over the broad X-ray band ~0.5 - 
-270 keV. The PSR J0205 +6449 X-ray spectrum above 2 keV 
has the same power-law shape (r ~ 1.03) as the middle-aged 
Vela pulsar, but the overall high-energy spectral shape, consid- 
ering also the Fermi y-ray spectrum, resembles more the spec- 
trum expected for a younger pulsar, i.e. no evidence for a thermal 
black-body component, and maximum luminosity at MeV ener- 
gies and not at GeV energies. 

The morphology of the double -pulse PSR J0205 +6449 light 
curve can be explained in a conventional outer-gap scenario for a 
rotating dipole in vacuum assuming low-altitude radio emission, 
similar to the case of the Crab pulsar when taking the Crab pre- 
cursor radio pulse as the counterpart of the single radio pulse de- 
tected for PSR J0205+6449. This can be verified in the "Atlas " 
of model y-ray light curves simulated by Watte rsetal. (2009). 
Howeve r, see also the alternative Atlas by Bai & Spitkovskv 
(2009a|), who point out an in consistency in the model calcu- 
lations by IWatters et al.l d2009l) and in earlier reports, affecting 
particularly profile shapes calculated for the two-pole caustic 
model. 

Furthermore, it should be realized that the sharp pulses in 
the high-energy profile of this young pulsar PSR J0205 +6449 
do not have radio counterparts like we see for the Crab, and 
we encourage a search for giant radio pulses in the phase in- 
ter vals of these high-energy pulses. Recent model calculations 
bylBai&Spitkovsky (2009b) using a force-free field instead of 
the vacuum dipole field show that alternative scenarios such as 
their annular-gap model are required to produce over a wide 
range of parameters two sharp high-energy pulses as exhibited 
by PSR J0205+6449. More extensive 3-D simulations including 
the physics of the production processes are required for more de- 
tailed comparisons with the spectral and timing characteristics. 
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